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Abstract: We present the NLO QCD corrections to pp — > 77 j production at hadron 
colliders. Our calculation includes contributions from the fragmentation of a hadronic jet 
into a highly energetic photon, and consequently allows the implementation of arbitrary 
infrared-safe photon isolation definitions. We compare different photon isolation criteria 
and perform a detailed study of the dependence of the 77J cross section on the photon 
isolation parameters. 
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1. Introduction 

Diphoton final states have played a crucial role in the recent discovery of a new boson at 
the LHC [1,2], and the fact that its partial decay width to photons seems to be enhanced 
as compared to the Standard Model (SM) Higgs boson keeps up the interest in this chan- 
nel. Diphotons are also important in many New Physics searches [3-6], in particular the 
search for extra spatial dimensions or cascade decays of heavy new particles. In particular, 
diphotons in combination with jets and missing energy occur in gauge mediated SUSY 
scenarios. 

In what concerns the SM processes at the LHC, the diphotons will most certainly be ac- 
companied by one or more high-p^ jets, which are often suppressed by a jet veto. However, 
the presence of an extra jet offers better control on the separation of backgrounds and sig- 
nals and more information on the interaction dynamics. Therefore a precise understanding 
of SM processes yielding diphotons in association with jets is mandatory. 

Diphoton production has been calculated at NLO some time ago [7], supplemented 
also by gluon initiated subprocesses beyond the leading order [8] and soft gluon resum- 
mation [9,10]. Recently, NNLO corrections to direct diphoton production also have be- 
come available [11]. Monte Carlo approaches to prompt photon production with matrix- 
element/parton-shower merging also have been studied [12-14]. In [12], the fragmentation 
contribution is generated by a combined QCD+QED shower, and photon isolation is done 



- 1 - 



within the democratic clustering approach [15-17]. Refs. [13,14] contain an NLO Monte 
Carlo study of diphoton production, where fixed cone-based photon isolation is also possi- 
ble. 

NLO calculations of diphoton-plus-jet production [18], photon-plus-two-jet produc- 
tion in QCD production [19] and in vector boson fusion [20], and diphoton production at 
NNLO [11] were up to now only carried out based on a smooth-cone isolation and did not 
admit the implementation of alternative isolation criteria. 

In this article, we present an NLO calculation of diphoton production in association 
with one jet which also contains a component from the fragmentation of QCD partons 
into photons, thereby allowing to study the effect of different photon isolation criteria. In 
addition, we provide a public code to compute 77+jet final states at NLO, where the virtual 
amplitude has been pre-generated with the automated one-loop program GoSam [21]. 
Section 2 contains a general discussion of photon isolation criteria and fragmentation. In 
Section 3 technical details of the calculation are discussed, while the results are presented 
in Section 4, before we conclude in Section 5. 

2. Photon isolation criteria and photon fragmentation 

Photons in hadronic final states can have multiple origins. Besides the radiation of a hard 
photon off a quark involved in the hard interaction (sometimes called 'prompt' photon), 
photons can also be emitted during the hadronization phase of a hadronic jet or can be the 
result of electromagnetic decays of unstable hadrons (these are sometimes called 'secondary 
photons'). Especially the decay 7r° — > 77 of pions at high transverse momentum can mimic 
the signature of a single photon if the two photons are too closely collimated to be resolved 
individually. Only the production of prompt photons can be computed within perturbation 
theory from first principles, while the production of photons in hadronization and hadron 
decays can only be modeled, thereby introducing dependence on ad-hoc parameters. 

Photons resulting from the hard interaction are usually well-separated from all hadronic 
jets produced in the event, while photons from hadronization and decay processes will 
always be inside hadronic jets. To disentangle prompt from secondary photons, one ap- 
plies isolation requirements, which limit the hadronic activity around a photon candidate, 
thereby defining an "isolated photon" . A veto on all hadronic activity around the photon 
direction would result in a suppression of soft gluon radiation in part of the final state 
phase space, thereby violating infrared safety of the observables. Consequently, all photon 
isolation prescriptions must admit some residual amount of hadronic activity around the 
photon direction. 

By admitting some hadronic activity around a photon, one includes final state con- 
figurations with a final state quark radiating a highly energetic collinear photon. These 
configurations contain a collinear singularity, related to small invariant masses of the quark- 
photon system. Mass factorization in QCD relates this singularity to a redefinition of the 
quark-to-photon fragmentation function, which describes the production of a photon in- 
side a hadronic jet. Like parton distributions in the proton, these fragmentation functions 
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are non-perturbative objects that have to be determined from experimental observations. 
Their dependence on the resolution scale is governed by evolution equations. 

The fragmentation contribution intertwines the production of prompt and secondary 
photons. An observable with final state photons defined through isolation criteria will 
typically receive some contribution from photon fragmentation, with the isolation crite- 
rion aiming to minimize this contribution. Several isolation criteria were proposed in the 
literature and applied in experimental studies of single photon and photon pair production. 

The cone-based isolation is most commonly used especially at hadron collider exper- 
iments. In this procedure, the photon candidate is identified (prior to the jet clustering) 
from its electromagnetic signature, and its momentum direction (described by transverse 
energy Etj, rapidity 7? 7 and polar angle 7 ) is determined. Around this momentum di- 
rection, a cone of radius R (typically chosen around 0.4) in rapidity 77 and polar angle (j) is 
defined. Inside this cone, the hadronic transverse energy £/r,cone is measured. The photon 
is called isolated if £r,cone is below a certain threshold, defined either in absolute terms, or 
as fraction of £r, 7 (typically 0.1 or below). The latter criterion then means that a photon 
candidate is considered as isolated if in a cone defined by 



the amount of hadronic energy does not exceed a pre-defined fraction of the photon energy: 



An alternative to the cone-based isolation is the democratic clustering procedure [15]. 
In this procedure, a photon candidate is treated like any hadron in the jet recombination, 
which can be performed using any (cone-based or clustering based) infrared-safe jet algo- 
rithm [22]. After the jet recombination, jets containing photons are labeled as photon jets, 
which are identified as isolated photons if the photon energy in the jet exceeds a predefined 
fraction z cut of the jet energy. Typical values of z cut are 0.9 or above. 

Both cone-based isolation and democratic clustering admit some fraction of events 
involving collinear quark-photon systems. The theoretical predictions for cross sections 
defined in these procedures must therefore take account of photon fragmentation contribu- 
tions. 

Finally, the smooth cone isolation criterion [23] varies the threshold on the hadronic 
energy inside the isolation cone with the radial distance from the photon. It is described 
by the cone size R, a weight factor n and an isolation parameter e. With this criterion, 
one considers smaller cones of radius r 7 inside the i?-cone and calls the photon isolated if 
the energy in any sub-cone does not exceed 



In theoretical studies, typical values used for the smooth cone isolation parameters are 
R = 0.4, e = 0.5, n = 1. By construction, the smooth cone isolation does not admit 
any hard collinear quark-photon configurations, thereby allowing a full separation of direct 
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and secondary photon production, and consequently eliminating the need for a photon 
fragmentation contribution in the theoretical prescription. Despite this advantage, the 
smooth cone isolation was up to now used in experimental studies of isolated photons only 
in a discretized approximation [24]. Owing to finite detector resolution, an implementation 
will only be possible to some minimal value of r 7 , thereby leaving potentially a residual 
collinear contribution. 

Perturbative calculations of isolated photon production must take proper account of 
the isolation criterion used to define the observable. Besides the usual higher-order QCD 
effects from unresolved partonic radiation, these calculations must also take account of the 
quark-photon collinear singularity appearing in the photon isolation procedure. In some 
specific observables (inclusive photon production in e + e~-annihilation [25, 26] and deep 
inelastic scattering [27-30] as well as photon-plus-one-jet production in e + e~ [15,16,31,32] 
and photon-plus- no-jet production in deep inelastic scattering [30]), this singularity appears 
already at the leading order [15,27]. These observables are consequently most sensitive to 
the photon fragmentation function and can thus be used for its determination [17,33]. 

Concerning the counting of perturbative orders, some ambiguity arises when photon 
fragmentation contributions are present. The photon fragmentation function itself is 0(a). 
Its evolution equation differs from the evolution equations for the fragmentation functions 
of hadrons [34] by a direct term, which results in logarithms of the resolution scale not 
being suppressed by one order of the strong coupling constant (as in the case of hadron 
fragmentation) . Motivated by this difference, it has been argued that the photon fragmen- 
tation function should be assigned an inverse power of the strong coupling constant [35-37], 
thereby shifting the relative order of direct and fragmentation contributions. Viewed by 
mass factorization counter terms (and infrared finiteness of the observables), the photon 
fragmentation function does not require an inverse power of the strong coupling constant. 
In calculations of electroweak corrections to jet observables [38-41] where similar photon 
isolation issues appear, the photon fragmentation function is consistently taken as O(a). 

Single photon and diphoton production at hadron colliders receive contributions from 
collinear quark-photon splitting only at next-to-leading order in perturbation theory. Among 
the existing calculations of NLO corrections to photon pair production [7-9] and photon- 
plus-jet production [42,43], the ones in [7,42,43] have implemented both cone-based and 
smooth cone isolation and assume the power counting of the photon fragmentation function 
to contain an inverse power of the strong coupling constant, which means that a multitude 
of fragmentation subprocesses have to be (and have been) included at NLO. In general, 
these subprocesses are however of rather minor impact on the total result after isolation 
cuts, owing to the overall smallness of fragmentation contributions in these observables. 
For the purpose of our calculations, we will therefore use a power counting of the photon 
fragmentation function as O(a). 

3. NLO corrections to pp — y 77+jet final states 
3.1 Structure of the calculation 

The calculation of the NLO corrections to pp —y 77+jet requires the combination of the full 
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QCD corrections with counter terms regularizing infrared QED singularities. In a first step 
we produced a code that is able to calculate the QCD corrections. For generation of the 
tree level and real emission matrix elements we use MadGraph [44,45], the regularisation of 
infrared QCD singularities is handled by MadDipole [46,47], which makes use of the dipole 
formalism as developed in [48]. For integration over the phase space we used MadEvent [49]. 
The routines for generating histograms and distributions originate from the MadAnalysis 
package (see http://madgraph.hep.uiuc.edu). 

The generation of the various pieces and their combination for the phase space inte- 
gration has been done in a fully automated way. 

3.2 Calculation of the virtual corrections 

The virtual corrections have been calculated with the automated one-loop amplitude gen- 
erator GoSam [21]. The program package GoSam starts from an input card edited by the 
user and generates the diagrams and the corresponding expressions for the loop amplitudes 
in an automated way, using QGRAF [50], FORM [51,52] supplemented with Spinney [53] for 
the spinor algebra, and haggies [54] for the automated code generation. It combines 
unitarity-inspired integrand reduction techniques [55-59] with traditional tensor reduction 
methods [60,61]. The rational part can be calculated algebraically within GoSam in an 
automated way. 

We generated two versions of the virtual contributions, with and without top loops 
and found the effects of virtual top quarks negligible. Therefore the results presented 
here have been obtained without the top contributions. Furthermore, we neglect one loop 
contributions with two initial state gluons, which are formally of higher order, as there is 
no corresponding tree level amplitude. Their contribution may potentially be enhanced by 
the large gluon luminosity [7,8]. Their inclusion is however beyond the scope of this paper. 

The remaining contributions can be reduced to the virtual corrections for the process 

qq^jjg- (3.1) 

The complete set of virtual corrections can be obtained from eq. (3.1) by crossing of the 
momenta and/or changing overall factors for different electromagnetic charges. In total 
there are 130 diagrams to be calculated up to pentagons. The virtual amplitude has 
been checked for gauge invariance by adding a momentum dependent part to the photon 
polarisation vectors. 

3.3 Calculation of real radiation contributions 

As explained in detail in Section 2 above, the QCD corrections for processes involving 
photons contain infrared singularities related to the collinear emission of the photon off 
a final-state QCD parton. These singularities are compensated by the mass factorization 
terms of the photon fragmentation functions. 

From the computational point of view this implies that the real emission matrix element 
contains collinear singularities which need to be regularized by some kind of subtraction 
terms. The corresponding integrated subtraction terms make this singularity apparent as 
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they develop an explicit pole term when integrated over the unresolved one-particle phase 
space in dimensional regularisation. This pole is then absorbed into the fragmentation 
functions. To regulate these singularities we again make use of the dipole formalism as 
developed in [62] and implemented in the QED extension of MadDipole [63]. This extension 
also offers the framework for a straightforward implementation of fragmentation functions. 
We refer to [63] for further details. 

A collinear singularity between a photon and a quark can be regulated by a single 
dipole. Note that in principle the role of emitter and unresolved particle is interchanged. 
In a true QED calculation the quark would be the emitter and the photon the unresolved 
particle. However in this calculation the photon is a tagged particle whereas the additional 
jet can be unresolved. Nevertheless the dipole formalism for QED can be used as the 
subtraction terms are symmetric under the exchange of emitter and unresolved particle if 
both are massless. A small modification has been made in such a way that also photons are 
allowed as spectators. This ensures that one can always use final-final dipole configurations, 
thereby reducing the complexity of the calculation. 

Upon phase space integration, the pointwise cancellation of the infrared poles from the 
virtual amplitude with those from the real radiation part has been checked. For both QCD 
and QED subtraction terms we have checked the independence on the a-parameter [64], 
which restricts the phase space segments on which dipole subtraction is performed to the 
vicinity of the infrared singularities. To perform this check we had to extend the usage of 
the a-parameter to the QED subtraction terms for non-collinear safe observables. In the 
purely massless case and for final-final configurations this extension is straightforward. 

NLO corrections to 77+jet final states for the Frixione isolation criterion are free 
of collinear quark-photon contributions and do not depend on the photon fragmentation 
function. For this specific isolation criterion, NLO results for 77+jet production have been 
derived in [18]. We fully reproduce these results, thereby obtaining a strong check on the 
correctness of our implementation of the virtual and real matrix elements, as well as on 
the non-QED-type subtraction terms. 

4. Numerical results 

In this section we present phenomenological results for proton-proton collisions at y/s = 
8TeV. The results are divided into two categories: the one-jet inclusive case, pp — > 
77+jet+X, and the one-jet exclusive case pp — > 77+jet. The inclusive sample is de- 
fined by requiring that at least one jet in the event passes the selection cuts introduced 
below, while the exclusive sample admits only events containing exactly one jet within the 
selection cuts. 

4.1 Input parameters and kinematic cuts 

For the jet clustering we used an anti-fcy algorithm [65] with a cone size of Rj = 0.4 
provided by the Fast Jet package [66,67]. We used an NLO parton distribution set from 
NNPDF2.3 [68], where the values for a s at leading order and next-to- leading order are 
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given by 

a s {M z ) = 0.119 , 

and the running is calculated at one loop for the LO results and at two loops for the 
NLO results. As we neglect initial state 6-quarks, we use the Nf = 4 version of the 
parton distribution set. For the photon fragmentation functions, we take set II of the 
parametrisations of Ref . [37] . 

The following kinematic cuts have been applied: pj? > 40 GeV, pj. > 20, |r/ 7 , rf\ < 2.5, 
Ryj > 0.4, i? 7i7 > 0.8 and 100 GeV < m 77 < 140 GeV. The intention of the latter cut is 
to focus on a region around the Higgs resonance. 

For the photon isolation, we compare the Frixione isolation criterion with the fixed 
cone criterion for several values of the photon energy fraction z c in the cone, where 

I rfh&d I 
I^T.cond 

Wt + PT,cone I 

such that in the collinear limit, z c is related to e c in eq. (2.1) by z c = j^-- For the Frixione 
isolation criterion (see eq. (2.2)), our default values are R = 0.4, n = 1 and e = 0.5. For 
the cone-based isolation, the default cone size is R = 0.4, while several different values of 
the hadronic energy fraction z c inside the cone will be used. 

4.2 Scale dependence and sensitivity to the isolation parameters 

The truncation of the perturbative expansion of a collider observable leads to a dependence 
on scale parameters that were introduced in the renormalization and mass factorization. 
The residual dependence on these parameters is often used to quantify the uncertainty on 
the calculation from missing higher order terms in the perturbative series. Apart from 
the dependence on renormalization and initial state factorisation scales /U r ,/i/, the cross 
section for the production of prompt photons also depends on the fragmentation scale hf, 
as explained in section 2 above. To study the scale dependence of our NLO results, we set 
H r = Hf = Hf and choose /Xq = \ (m 77 + J2j Pt j) f° r our cen t ra l scale. The scales are then 
varied by \x = x fio with 0.5 < x < 2. 

In Figure 1, we display the scale dependence of the exclusive 77+jet cross section. 
For a cone-based isolation criterion, we observe a clear reduction of the scale dependence 
at next-to-leading order, while this reduction is less pronounced for the Frixione isolation 
criterion. This qualitative difference can be attributed to the occurrence of almost collinear 
quark-photon configurations at NLO. The typical scale of these configurations is the in- 
variant mass of the quark-photon system, which can be substantially lower than ^o- I n the 
case of a cone-based isolation, these contributions combine with the fragmentation con- 
tribution, which compensates their scale dependence. For the Frixione isolation criterion, 
this compensation does not occur, thereby resulting in a larger scale-dependence. 

In the inclusive case, Fig. 2, no reduction of the scale uncertainty at NLO is visible. The 
reason for this is the fact that the cross section in this case is dominated by the 77+2 jets 
real radiation part, which shows a leading order scale dependence. A similar behaviour has 
been observed for example in ZZ+jet production with and without second jet veto [69]. 
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The largeness of the NLO corrections has already has been observed for Frixione isolation 
in Ref. [18]. The reasons are the appearance of new partonic subprocesses at NLO and the 
enlarged final state phase space at this order. 
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Figure 1: Behaviour of the exclusive 77+jet cross sections with different isolation prescriptions 
under scale variations, /j, = x /^o, 0.5 < x < 2, (1q = j (m^ 7 + X^Ptj)- 




Figure 2: Behaviour of the inclusive 77+jet+X cross sections with different isolation prescriptions 
under scale variations, [i = x fiQ, 0.5 < x < 2, /Iq — i + ^2jPr j)- 

Due to photon isolation, the cross section for pp — > 77+jet+X is not strictly an in- 
clusive quantity. The integration over the final state collinear variable z appearing in the 
fragmentation functions Dq(z, /if) is bounded from below by 1 — z c . Further, the presence 
of three different scales fj, r , fj,f, fj,p partially leads to a scale dependence which is different 
from what is known from pure QCD. 

Figures 3 and 4 show how the scale variation bands vary as a function of the isolation 
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parameters, for both the single-jet inclusive and the exclusive case. One observes that 
for both isolation criteria the inclusive case is dominated by the large scale dependence of 
the 77 + 2 jets part of the real radiation which has an uncompensated leading order scale 
dependence. 

3 
2.8 
2.6 
2.4 
2.2 

a. 

" 2 
o 

6" 

1.8 
1.6 
1.4 
1.2 
1 

0.02 0.04 0.06 0.08 0.1 0.12 0.14 0.16 0.18 0.2 

z c 

Figure 3: Dependence of the cross sections on the cone isolation parameter z c . The bands corre- 
spond to scale variations 0.5 < x < 2, with /i = x /io. 
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Figure 4: Dependence of the cross sections on the Frixione isolation parameter e. 

The small size of the scale variation band for z c < 0.14 in the case of exclusive cuts 
with cone isolation can be understood from the fact that in this range of z, the log(// j p) 
terms which are contained in the perturbative component of the fragmentation functions 
largely cancel the [ip dependence of the direct real radiation part. Only for larger z c values 
non-perturbative and beyond-leading-logarithmic effects start to become important. We 
validated this cancellation in detail on the NLO corrections to diphoton production, as 
implemented in the code diphox [7]. 

With Frixione isolation, the log(//i?) dependence is absent, therefore the behaviour 
under scale variations in Fig. 4 is qualitatively different. 
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4.3 Results for diphoton plus one jet production — exclusive case 

In the exclusive case, only moderate NLO corrections are observed in the kinematical 
distributions related to 77+jet final states. Figure 5 displays the photon pair invariant 
mass and the jet transverse momentum distributions. We observe that the NLO corrections 
are in general limited in magnitude, and slightly larger for the Frixone isolation than for 
the cone-based isolation. Inclusion of the NLO corrections amounts to an rescaling of the 
distributions that is constant in invariant mass, and slightly decreasing with jet transverse 
momentum. 




M„[GeV] PTje.,[ GeV ] 



Figure 5: (a) Photon invariant mass distribution, (b) transverse momentum distribution of the 
leading jet for the diphoton plus one jet exclusive cross section. 



In Figure 6, the transverse momentum distributions of the leading and the subleading 
photon are displayed. It can be seen that the leading-p^ photon tends to become softer at 
NLO, which is natural since the additional QCD radiation carries away momentum. 
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Figure 6: Transverse momentum distributions of the leading (i.e. larger-p^) and subleading photon 
for the diphoton plus one jet exclusive cross section. 

Figure 7 shows the differential distributions in the distance Rj^ = \J {r\i — r/T) 2 — (r£? — </> 7 ) 2 
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between the jet and the harder photon (71) respectively the softer photon (72), with a sepa- 
ration cut of Rj^y > 0.4. At leading order, the preferred kinematical configuration of the jet 
is back-to-back with the harder photon, and near to the softer photon. These kinematical 
correlations are weakened with the inclusion of NLO corrections, with an increase of events 
with smaller opening angle between the hard photon and the jet or larger opening angle 
between the soft photon and the jet. 




let 1,7 1 R |et1, r 2 

Figure 7: i?-separation between photon and jet in the r\ — <j) plane for the diphoton plus one jet 
exclusive cross section. Rj etini denotes the i?-separation between the jet and the harder photon, 
while Rj etl , l2 is the i?-separation between the jet and the softer photon. 
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4.4 Results for diphoton plus one jet production — inclusive case 

As already observed for the total cross section, NLO corrections are substantially larger for 
the inclusive cross section 77+jet+X as compared to the exclusive case. In the inclusive 
cross section, the substantial contribution from 77 + 2 jet final states results in larger 
corrections, and induces substantial modifications to some of the kinematical distributions. 

Figure 8 displays the inclusive distributions in photon pair invariant mass and leading 
jet transverse momentum. The magnitude of the corrections is larger than in the exclu- 
sive case, they remain constant for the invariant mass distribution and rise with the jet 
transverse momentum (as opposed to the decrease with jet transverse momentum in the 
exclusive case, Figure 5). Again, the corrections for the Frixione isolation criterion are 
slightly larger than for the fixed-cone isolation. 
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Figure 8: (a) Photon invariant mass distribution, (b) transverse momentum distribution of the 
leading jet for the diphoton plus one jet inclusive cross section. 

The photon transverse momentum distributions, Figure 9, display a similar behaviour 
as in the exclusive case, with the main effect from NLO corrections appearing in a softening 
of the leading photon distribution. The effect of the extra jet in the inclusive distribution 
is particularly pronounced in the Rjy distributions. 

Comparing exclusive (Fig. 7) and inclusive (Fig. 10) cases, one can see very clear 
differences. For example, in the first bins of the Rj^ distribution (separation between 
leading jet and harder photon), the inclusive case shows a shoulder due to the contributions 
from the second jet, which is vetoed in the exclusive case. Further, in the first bins of the 
Rj 12 distribution (separation between leading jet and softer photon), the K-factor is smaller 
than one in the exclusive case, while it is always larger than unity in the inclusive case. 
Note however that in the inclusive case, events where both jets fulfill the cuts are counted 
twice. Therefore it is somewhat misleading to directly compare K-factors between the 
inclusive and exclusive case. 
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P Tjl [GeV] P Tj2 [GeV] 

Figure 9: Transverse momentum distributions of the leading (i.e. larger-p^) and subleading photon 
for the diphoton plus one jet inclusive cross section. 




R |et1, 7 1 R |el 1,-, 2 

Figure 10: i?-separation between photon and jet in the rj — cf> plane for the diphoton plus one jet 
inclusive cross section. Rjet 1 ,-y 1 denotes the -/^-separation between the jet and the harder photon, 
while i?jet!,7 2 is the i?-separation between the jet and the softer photon. 
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5. Conclusions and outlook 

We have calculated the NLO QCD corrections to the production of two isolated photons in 
association with a jet at hadron colliders. Results for both the one-jet inclusive case as well 
as the case where exactly one jet passing the jet selection criteria have been presented. Our 
calculation also includes contributions from the fragmentation of a hadronic jet into a large- 
ly photon at order a 2 a 2 s , and therefore allows to compare different photon isolation criteria. 
Comparing the Frixione isolation criterion to standard cone isolation for several values of 
the hadronic energy allowed in the cone, we observe that the scale dependence stabilizes at 
NLO for the exclusive cross section in the standard cone isolation case, while with Frixione 
isolation the stabilisation is much less pronounced. This behaviour can be attributed to the 
fact that there are cancellations of logarithms stemming from the factorisation of collinear 
quark-photon splittings between the direct and the fragmentation contributions. For very 
strict isolation parameters, the results for standard cone and Frixione isolation are similar 
with regards to the scale dependence. The K-factors are in general larger with Frixione 
isolation than with cone isolation. 

In contrast to the exclusive cross section, the one-jet inclusive cross section does not 
show a stabilisation of the scale dependence, independent of the choice of the isolation 
criterion. This can be understood from the fact that the cross section in this case is 
dominated by the two-jet contribution of the NLO real radiation part, which shows a 
leading order scale dependence. 

The code which is underlying the calculation presented here will be made publicly 
available at http://gosam.hepforge.org. This will allow further dedicated studies of 
different isolation prescriptions and kinematic situations, and eventually the combination 
of the NLO calculation with a parton shower. 
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